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ABSTRACT 


An  ideal  turbomachine  tranters  energy  from  a  rotating  shaft  to  a  fluid  through  the 
creation  of  a  flow  field  that  is  unsteady  in  the  absolute  frame.  In  order  to  understand  and 
control  the  work  and  acoustics  of  this  process  the  basic  building-block  singularities  of 
potential  flow  are  developed  in  a  unified  manner  such  that  both  near-  and  far-field  effects 
appear  naturally  in  the  formulation.  It  is  shown  how  the  total  vector  force  applied  to  the 
fluid  by  a  discrete  rigid  blade  is  related  to  the  bound  and  shed  vortex  system  and  the  ac¬ 
celeration  of  the  displaced  mass.  Following  Lighthill,  the  volume  integral  of  the  vector 
moment  ofvorticity  is  found  to  describe  most  fundamentally  how  the  unsteady  fluid  mo¬ 
tions  in  a  turbomachine  are  related  to  surface  forces  that  do  work  on  the  fluid. 

Because  unsteady  flow  is  required  for  energy  transfer,  the  work  input  of  a  single 
blade  creates  simultaneously  in  the  fluid  an  acoustic  signal  that  propagates  to  the  far- 
field.  The  strength  of  this  signal  is  determined  by  the  time  derivative  of  the  instantaneous 
blade  force  component  in  the  direction  of  a  stationary  observer,  i.e.,  with  the  pattern  of  a 
rotating  dipole.  The  summation  of  these  tonal  signals  from  combinations  of  blades  and 
vanes  is  computed  for  free-field  and  cylindrical  duct  boundary  conditions.  Bessel  func¬ 
tions  are  found  to  characterize  both  cases,  their  order  being  determined  by  specific 
blade/vane  number  differences.  Rules  for  tonal  sound  control  by  acoustic  interference  or 
mode  decay,  such  as  those  of  Tyler  and  Sofrin,  are  examined.  Implications  for  broadband 
acoustics  are  discussed. 

The  ideal  unsteady  singularity  concepts  are  exteruled  to  turbomachines  having  cen¬ 
trifugal  and  real  fluid  effects.  It  is  concluded  that  this  unified  vortexidipole  approach 
provides  a  sound  basis  for  continued  turbomachinery  understanding  and  de\>elopments. 

ADMINISTRATIVE  INFORMATION 

This  report  was  prepared  by  the  research  staff  of  the  Propulsion  and  Auxiliary  Sys¬ 
tems  Department  (Code  27)  of  the  David  Taylor  Research  Center  (DTRC)  for  the  purpose 
of  organizing  the  theoretical  knowledge  relating  to  a  machinery  component  common  to 
several  department  projects.  It  was  supported  by  department  overhead  funds. 

INTRODUCTION 

The  mathematics  of  potential  flow  represent  a  simple,  yet  effective,  analysis  tool 
that  can  provide  an  accurate  initial  understanding  of  many  three-dimensional  time-vary¬ 
ing  flows.  Moreover,  this  technique  has  the  desirable  property  that  solutions  are  additive 
so  that  particular  phenomena  can  be  studied  separately  in  simple  form,  and  then  summed 
to  represent  complex  cases.  One  such  case  is  the  turbomachine,  which  is  made  up  of  a  set 
of  discrete  air-foil  shaped  blades  rotating  in  a  three-dimensional  flow  field.  From  the 
thermodynamic  stand  point  this  machine  receives  input  through  a  shaft  that  rotates  at 
(generally)  constant  speed  and  experiences  a  (generally)  constant  torque  and  thrust.  As  a 
result,  the  machine  adds,  or  subtracts,  work  from  the  moving  fluid,  and  may  also  give 
propulsive  thrust  to  a  vehicle.  Using  the  building  block  property  of  potential  flow  theory, 
two-dimensional,  three-dimensional,  and  unsteady  flows  can  be  examined  separately  and 
in  combination.  For  turbomachines  it  is  especially  revealing  to  simulate  the  motion  of  the 
blades,  and  the  resulting  time-varying  velocity  fields,  as  seen  by  a  fixed  observer  in  order 
to  account  mathematically  for  the  work  done  on  the  fluid.  But,  beyond  the  work,  and  as  a 
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usually  unwanted  side  effect,  the  blade  motion  also  initiates  acoustic  signals  in  the  fluid, 
which  propagate  at  the  speed  of  sound  over  great  distances.  Since  most  acoustic  analyses 
are  based  on  potential  flow  theory,  it  is  the  goal  of  this  study  to  develop  a  potential  flow 
description  of  the  fundamental  machine  work  process  that  leads  directly  to  a  unified  tech¬ 
nique  for  characterizing  the  sound  field  associated  with  the  work. 

The  concepts  developed  herein  are  an  extension  of  an  earlier  study  by  this  author. 
That  study  demonstrated  the  advantage  of  using  unsteady  potential  flow  vortex  fields  to 
describe  the  work  and  acoustic  signature  of  an  axial  flow  rotor  in  a  three-dimensional  free 
field.  Since  that  initial  investigation,  it  has  been  found  useful  to  explore  further  the 
theory,  to  consider  centrifugal  effects,  and  to  include  the  details  of  acoustic  interference 
of  rotor-stator  combinations  in  open  and  ducted  environments.  In  applying  the  theory,  this 
smdy  builds  on  the  fluid  dynamic  principles  presented  by  Lamb,  since  that  is  the  classic 
treatment  of  potential  flow  and  basic  acoustics.  Another  very  helpful,  and  contemporary, 
source  of  background  is  the  work  of  Lighthill,  both  for  potential  flow  and  acoustics.  Oth¬ 
er  singular  contributions  are  also  used,  such  as  Preston  for  the  moving  vortex -work 
analysis,  Gutin  for  the  original  acoustic  interference  theory,  lyier  and  Sofrin  for  the  rotor- 
stator  acoustics,  and  Goldstein,  Cumpsty  and  Blake  for  their  comprehensive  summaries 
of  the  Aero-Hydro  Acoustics  of  turbomachines. 

It  is  therefore  the  purpose  of  this  report  to  present  the  essential  aspects  of  unsteady 
potential  flow  vortex  systems  so  as  to  provide  a  complete  and  consistent  description  of 
the  work  and  acoustic  behavior  of  low  speed  turbomachines.  The  condition  of  “low 
speed”  is  used  both  to  eliminate  the  need  to  consider  significant  fluid  compressibility,  and 
to  require  that  acoustic  signals  generated  over  any  given  time  interval  will  have  propa¬ 
gated  distances  much  greater  than  a  blade  dimension  or  movement  in  that  same  time.  In 
addition,  all  effects  are  considered  to  occur  without  fluid  viscosity,  without  turbulence, 
and  with  rigid  stationary  and  moving  surfaces.  The  following  pages  describe  bow  the 
work-related  parts  of  turbomachine  time-varying  velocity  fiel^  are  viewed  best  as  ex¬ 
panding  three-dimensional  vortex  loops,  or  equivalently,  to  injected  dipoles;  how  the 
important  features  of  displaced  and  added  mass  arise  for  an  oscillating  rigid  sphere  di¬ 
pole;  how  multisource  acoustic  interference  develops  from  multiple  blades  in  a  ffee-field 
and  in  a  duct  environment;  and  bow  the  basic  vortex/dipole  concept  can  give  insight  into 
even  more  complex  turbomachine  behavior. 

BLADE  VORTEX  LOOPS,  DIPOLES  AND  WORK 

THE  UNSTEADY  VELOCITY  POTENTIAL  AND  ENTHALPY 

In  order  to  apply  the  mathematics  of  potential  flow  it  is  necessary  to  assume  that  the 
fluid  have  no  viscosity  and  be  in  irrotational  motion,  i.e.,  =  0  and  the  vorticity,  ^  =  0 

throughout  the  flow  field.  These  two  restrictions  are  acceptable  for  many  flows  and  are 
often  related  since  viscosity  generates  rotation  in  an  otherwise  irrotational  field.  When 
these  conditions  are  met,  and  the  fluid  is  essentially  incompressible,  it  is  possible  to  ex¬ 
press  the  velocity  anywhere  in  the  flow  field  in  terms  of  the  gradient  of  a  scalar.  Thus  the 
equation  for  the  conservation  of  mass  in  the  fluid  becomes  Laplace’s  equation,  solutions 
to  which  require  only  boundary  conditions.  Hence, 

0  =  0(j:,  t)  and  V  =  V0,  then  V  •  V  =  0  =  V^<f>  (1) 
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Therefore  <p  is  termed  the  “velocity  potential”  of  a  flow  field,  with  the  convenient 
property  that  multiple  potentials  can  be  summed  to  represent  more  and  more  complex 
flows.  One  of  the  building  block  potentials  that  describes  circular  motions  is  that  of  the 
vortex  at  the  origin,  which,  for  a  two  dimensional  field  is  simply; 


(2) 


Here  d  is  an  angle  measured  from  some  reference  axis,  and  F  is  the  strength  of  the  vortex. 
Applying  Eq.  (1),  the  two-dimensional  velocity  components  are: 


Vr  =——  =  0  and  vg  =  — — 
dr  r  86 


Thus  the  radial  velocity  is  zero  and  the  circular  velocity  varies  inversely  with  the 
radius,  being  singular,  i.e.  infinite,  at  the  origin.  It  can  be  seen  that  a  line  integration  of 
the  velocity  around  any  closed  loop  that  includes  the  singularity  gives  the  value  F  .  Be¬ 
cause  of  this,  the  velocity  potential  of  the  vortex  can  be  seen  to  have  the  rather  unique 
property  that  it  increases  by  this  value  for  any  2ji  change  in  angle,  a  behavior  most  im¬ 
portant  for  simulating  turbomachine  force  and  work  effects. 

In  addition  to  the  conservation  of  mass,  the  conservation  of  linear  momentum  must 


also  be  satisfied  for  any  flow  system,  i.e.. 


8V  1  ?:/ vA 


.f  X  V  +  -^V^V 

Qo 


For  potential  flow  of  an  incompressible  fluid  this  equation  can  be  simplified  to  the  un¬ 
steady  Bernoulli  form. 


Therefore  a  gradient  in  fluid  total  enthalpy,  which  is  the  purpose  of  a  mrbomachine, 
is  associated  directly  with  a  gradient  in  the  time  variation  of  the  velocity  potential.  If  the 
vortex  of  Eq.  (2)  is  steady  in  time,  then  the  total  enthalpy  is  uniform  throughout  the  flow 
field,  giving; 


A  SINGLE  MOVING  2-D  VORTEX 

Usually,  the  static  pressure  at  r  =  oo  is  fixed,  at  say  zero,  then  the  pressure  at  the 
origin  will  be  -  oo  .  if  then  the  vortex  is  allowed  to  translate,  say  in  the  -x  direction  with 
velocity  Uj  ,  the  behavior  of  total  enthalpy  with  respect  to  the  absolute  frame  is; 
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d<b  T  J  y  \  T  Ut  ^ 

^  =  — tan-M  — ^  -i-syndrei  =  -  UjUrei  =  -Ht  (7) 

dt  2jt  \  A'  +  Uft  J  2n  Tret 

where  0,,/  and  r^.,/  are  measured  from  the  center  of  the  moving  vonex.  As  shown  in  Fig. 
1  there  is  now  a  change  in  the  fluid  enthalpy  when  viewed  in  the  absolute  frame,  as  one 
moves  across  the  x  axis. 


ng.i.  Contours  of  constant  total  enthalpy  for  a  moving  vortex. 

For  a  simple  two-dimensional  turbomachine  there  is  also  flow  across  the  vortex,  which 
can  be  described  by  additional  velocity  potentials  in  both  the  positive  x  and  y  directions: 

<pF=Uj,x+Uyy  (8) 

Looking  at  the  vortex  center,  in  a  frame  of  reference  attached  to  the  vortex,  there  is  thus  a 
relative  velocity. 

^rel  =  KUx  +  Uj)  +jUy  (9) 

If  one  computes  the  change  in  the  convection  of  fluid  momentum  across  planes  nor¬ 
mal  to  the  direction  of  the  relative  velocity  and  at  equal  distances  in  front  of  and  behind 
the  vortex  position,  the  result  is  the  force  applied  to  the  fluid  by  the  vortex  as  indicated  in 
Fig.  2. 

^^rel  ^  j  j  \dA  =  QoVrelh  ^  ^  Jdn  (10) 
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Fig.  2.  Boundaries  of  integration  for  computing  force  and  work 
produced  by  a  moving  vortex. 
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Examining  the  total  fluid  field,  it  is  seen  that  there  is  convection  of  mass  perpendic¬ 
ular  to  the  direction  of  vortex  motion  by  Uy.  If  we  compute  the  increase  in  energy  of  the 
mass  flow  between  any  two  planes  on  either  side  of  the  x  axis,  the  result  is,  using  Eq.  (5); 


J  QoUyft{HT2  -HT^)dx  =  QoUyh  j 


00 


00 


(15) 


^/Af  =  QoUyhTUT  =  Q^relh 


Ut  =  F,Ut 


(16) 


Hence  the  unsteady  vortex  mathematics  is  consistent  in  that  the  power  introduced  into  the 
fluid  by  the  moving  vortex  is  represented  exactly  by  the  rise  in  total  enthalpy  of  the  fluid 
convected  across  the  path  of  the  vortex.  Thus  the  solution  for  the  moving  single  vortex 
simulates  the  effect  on  the  fluid  of  a  moving  airfoil  at  some  angle  of  attack. 

The  global  energy  change,  after  some  time,  will  appear  in  the  wake  of  the  vortex 
around  a  line  defmed  by  the  relative  velocity.  The  reason  for  this  wake  development  is 
associated  with  the  very  important  fact  that  in  this  potential  flow  there  is  also  a  require¬ 
ment  to  conserve  vorticity.  Therefore  the  instant  that  the  work  vortex  is  initiated  a  second 
vortex  of  equal  strength,  but  of  opposite  direction,  is  created  simultaneously.  This  second 
vortex  is  convected  with  the  flow,  but  does  not  contribute  to  any  further  work  input  since 
it  has  no  force  on  it.  However,  as  shown  by  Wu,  the  velocity  potential  of  a  vortex  pair  is 
equivalent  to  a  line  of  uniform  volume  dipoles  stretched  between.  Thus  the  time  behavior 
of  this  line  dipole  accounts  for  the  energy  change  in  the  wake.  Figure  3  illustrates  the  2-D 
streamlines  for  this  situation. 


/ 
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A  ROW  OF  MOVING  2-D  VORTICES 

The  previous  analysis  suggested  how  a  single  vortex  interacts  with  a  fluid  to  ex¬ 
change  energy  in  a  two-dimensional  flow.  This  approach  has  been  generalized  to  an 
infinite  row  of  two-dimensional  vortices,  and  to  a  circular  arrangement  of  vortices,  by 
Preston.  These  two-dimensional  geometries  simulate  the  action  of  axial  flow  and  centrif¬ 
ugal  flow  mrbomachines  respectively.  Taking  for  example  the  infinite  row  of  vortices  of 
strength  F  and  separation  a,  Preston  adapts  a  steady  solution  for  the  velocity  potential 
from  Lamb. 

<b= — tan*^  (cot—  tanh  — ]  (17) 

^  2jc  \  a  a  ) 


If  this  vortex  row  now  moves  to  the  -x  direction  at  velocity  Ut  ,  then  for  some  Xo  in  the 
absolute  frame,  x  =  Xo  +  Ut  t,  and  the  potential  becomes  a  function  time.  Taking  the 
space  or  time  derivatives,  there  results  for  any  point  in  the  flow  field,  with  a  superim¬ 
posed  flow: 


dXo 


d0 

dy 


=  v=  t/v  +  - 


2a 


sinh  {2jty/a) 

cosh  i2jcy/a)  -  cos(2jtx/a) 
sin(2jtXo/a) 

cosh  (2ny/a)  -  cos(2jcc/a) 


(18a) 


(18b) 


Ai^.  =  -A(Hr)  =  —  sinCWa) 

dt  2a  cosh  (Tsvyja  -  cos{2jtXo  +  U  jt)/  a 


(18c) 


Looking  on  either  side  of  the  x  axis  gives  the  following  asymptotic  values  for  the  total 
enthalpy: 


H[(y<a)  = 


TUt 

2a 


(19a) 


Ht  (y<-a)  = 


~rur 

2a 


=  Uju. 


(19b) 


Therefore  the  change  in  enthalpy  of  a  fluid  particle  as  it  moves  across  the  row  of  translat¬ 
ing  vortices  is  Uj  F/a  or  Ut  ( 2u+),  If  the  effect  of  the  remote  starting  vortex  row  is 
included,  as  in  Fig.  4,  the  enthalpy  changes  are  the  same,  however,  the  velocity  change  is 
zero  for  y  <-  a  and  20^  for  y  >  a.  Although  illustrated  here  for  a  two-dimensional  carte¬ 
sian  geometry,  the  result  for  the  circular  arrangement  of  vortices  is  similar.  It  may  be 
noted  that  this  analysis  does  not  define  how  the  vortex  is  generated,  only  that  it  be  asso¬ 
ciated  with  a  force  on  the  fluid  having  a  component  in  the  direction  of  the  vortex  motion. 
For  the  potential  flow  description  of  a  turbomachine  it  is  asstuned  that  each  blade  is  asso¬ 
ciated  with  a  vortex,  or  a  group  of  vortices  of  variable  strength,  as  required  to  match  the 
blade  surface  boundary  condition.  If  the  vortices  do  not  move  they  represent  a  cascade  of 
stationary  turning  blades  where  the  total  pressure  is  constant,  but  the  static  pressure  may 
change  according  to  the  Bernoulli  relation,  Eq.  (6). 
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Fig.  4.  Bound  and  starting  vortex  arrays  for  a  linear  blade  row. 

VORTEX  LOOPS  AND  DIPOLES  IN  3-D 

For  three-dimensioiial  axial  flow  turbomachinery  the  vortex  velocity  potential  pro¬ 
duces  similar  effects;  however  for  this  case  it  is  necessary  to  examine  more  closely  the 
connection  between  the  blade  vortex  and  the  starting  vortex.  In  three-dimensions  the 
blade  has  only  a  finite  span,  so  the  vortex  behavior  at  the  root  and  the  tip  ends  must  be 
accounted  for.  This  is  accomplished  by  introducing  two  “shed”  vortex  lines  from  the  ends 
of  a  uniformly  loaded  blade,  each  having  the  same  strength  and  continuity  of  rotation  as 
the  blade  vortex.  These  shed  vortex  lines  are  convected  with  the  flow  through  the  ma¬ 
chine  and  eventually  connect  with  the  starting  vortex.  Thus  in  three-dimensions  the  basic 
turbomachine  vortex  system  becomes  1 )  a  line  vortex  along  the  blade  span,  2)  two  line 
vortices  convecting  downstream  from  the  hub  and  tip,  and  3)  a  starting  vortex  in  the  far 
downstream  region.  As  shown  in  Fig.  5,  this  complete  vortex  system,  trailing  in  a  spiral 
path  from  the  blade  disc,  will  describe  all  the  pressure,  velocity  and  enthalpy  changes  in¬ 
troduced  by  the  blade.  The  essential  features  of  the  vortex  loop  trailing  from  each  blade 
are  constant  strength,  and,  area  increasing  at  a  constant  rate  as  the  blade  rotates.  Consid¬ 
ering  the  situation  in  the  vicinity  of  the  blade  row,  it  is  the  growth  of  this  vca'tex  loop  that 
is  critical  to  the  work.  In  the  downstream  flow,  the  mathematical  combination  of  all  the 
spiralling  vortex  loops  from  many  blades  creates  in  the  fluid  a  region  of  increased  enthal¬ 
py  that  has  been  acted  on  by  the  blade  row,  and  outside  of  this  affected  region  there  are 
no  changes. 
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Rg.  5.  Three  dimensional  helical  vortex  loop  created  by  a  rotating 
blade  in  axial  flow. 


TXiming  to  the  fluid  Held  near  each  blade,  it  can  be  seen  that  the  vortex  loop  com¬ 
posed  of  the  blade  vortex  and  the  tip  and  root  shed  vortices  is  increasing  in  area  at  the 
rate: 

=  (2C 

dt 

Now,  Lamb  shows  that  a  vortex  loop  of  area  A  and  strength  F  in  an  undisturbed 
flow  is  equivalent  to  the  impulse  that  was  necessary  to  start  that  loop  in  motion. 

^VL  =  Qt^AvjJlA.  =  [  Fyidt  (212 


A  steady  increase  of  vortex  loop  area  is  then  equivalent  to  a  steady  force  applied  to 
the  fluid  such  that 

Fvl  =  Q(X  X  ^relh  (21b) 

which  is  the  three-dimensional  equivtdent  of  Eq.  (12).  Furthermore,  the  steady  velocity 
potential  at  any  remote  point  due  to  an  incremental  area  of  vortex  loop  at  the  origin  is 
given  by  Lamb  as 

,  r  ♦  ♦  AAvl 

<PvL  =  —  ■  /la— v-  (22) 

4jr  fo 
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where  rig  is  the  vector  from  the  origin  to  the  field  point  and  is  the  distance.  The  form 
of  (22)  is  such  that  it  is  equal  to  the  velocity  potential  of  a  steady  volume  dipole  located  at 
the  area  centroid,  and  having  a  direction  normal  to  the  plane  of  the  vortex  loop  area  as 
illustrated  in  Fig.  6,  and  given  by  (23). 


Fig.  6.  Coordinate  definitions  for  the  velocity  potential  of  a 
three-dimensional  vortex  loop. 


4>DirD,P) 


M>y^  * 


(23) 


Thus  the  incremental  dipole  strength  is  equivalent  to  FAA,  and  it  can  be  seen  that 
continuous  expansion  of  vortex  loop  area  is  equal  to  continuous  injection  of  dipoles  into 
the  flow  field  from  the  position  of  each  blade.  In  fact,  the  DNA-like  spiralling  double  he¬ 
lix  of  the  vortex  pair  shed  into  the  wake  of  each  blade  can  be  replaced  mathematically 
with  a  single  hehcal  line  dipole  centered  between  the  vortex  pair.  From  an  equivalence 
standpoint  it  is  possible  to  say  all  of  the  following: 


dly 


Fy  =  —  =  QgTxVreih  =  QgTX 

at 


^  dAvL 


~  Qo  j  —  Qo^rel  ~Z 
dt  ds 


DIPOLE  COMPONENTS  FROM  GREEN’S  THEOREM 

In  addition  to  the  work  done  by  a  blade  as  it  moves,  the  fact  that  each  blade  has  a 
finite  volume  means  that  fluid  must  be  displaced  continuously  from  the  front  and  pulled 
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into  the  back  as  the  blade  moves.  The  potential  flow  solution  for  the  case  of  a  rigid  body 
moving  through  a  fluid  is  presented  fully  by  Lamb.  However,  recently,  LighthiU  has  pres¬ 
ented  a  more  physically  significant  approach  to  this  situation.  Following  LighthiU,  the 
Green’s  surface  integral  theorem  allows  a  velocity  potential  at  any  position  in  a  flow  field 
to  be  computed  from  two  integrals  of  the  originating  potential  over  the  surface  of  the 
body  creating  the  flow.  When  the  point  p  is  geometrically  remote  from  the  body; 


Evaluation  of  these  integrals  for  the  case  of  a  general  body  shape  shows  that 

<tm,  ds  =  Gi=  I/Qo 

1 =  = 

where  I  is  the  impulse  applied  to  the  fluid  external  to  the  body  to  establish  that  flow 
field,  and  A  is  the  volume  of  fluid  displaced  by  the  body.  In  addition  to  the  above  effects 
of  rigid  body  motion.  LighthiU  shows  that  if  there  is  a  circulation,  F ,  associated  with  the 
body,  as  is  the  case  for  the  lifting  airfoil,  then  there  must  be  a  three-dimensional  vortex 
system  in  the  region,  and  that  the  lift  force  in  particular  is  given  by  the  time  rate  of 
change  of  the  volume  integral  of  the  vector  moment  of  vorticity  over  the  entire  region: 

=  “  1 1  X  f )  (27a) 


(26a) 

(26b) 


Rg.7.  Coordinate  system  for  evaluation  of  vorddty  moment  volume  integral. 
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Here  the  integral  must  be  taken  over  all  space  so  as  to  include  the  whole  vortex  loop.  Of 
course  in  potential  flow  the  bulk  vorticity  is  zero,  but  in  singular  regions,  where  the  vor¬ 
tex  lines  exist,  one  can  say: 


^dA  =  F/Zv 


(27b) 


so  that  the  integral  of  (27a)  gives  the  force  of  (24).  Now,  as  Lighthill  points  out,  it  just  so 
happens  that  the  impulse  term  of  (26a)  can  be  related,  through  a  theorem  by  Kelvin,  to  a 
system  of  vortex  loops  placed  in  the  surface  of  the  body  to  account  for  the  tangential  slip 
velocity.  If  then  these  vortex  loops  are  used  to  describe  the  impulse  term,  Lighthill  ob¬ 
tains  the  simple,  but  acoustically  useful,  result  for  the  total  force  to  move  a  blade  body 
through  a  fluid. 


It  should  be  noted  that  from  a  large  distance  all  of  the  surface  shear  vortex  loops  can  be 
combined  to  appear  as  one  dipole.  Interestingly,  for  harmonic  motion,  the  only  possible 
net  power  input  to  a  non-lifting  body  is  through  the  surface  shear  term,  so  it  will  be  illus¬ 
trative  to  study  later  in  more  detail  the  oscillating  motion  of  dipoles.  Before  examining 
periodic  dipoles  however,  it  is  important  to  discuss  the  vortex/dipole  system  associated 
with  centrifugal  flow  turbomachines. 


VORTEX  FORMATION  IN  CENTRIFUGAL  ROTORS 

The  previous  models  of  unsteady  potential  flow  were  developed  implicitly  for  linear 
or  axial  flow  machines  where  there  is  little  fluid  motion  in  the  radial  direction.  Although 
it  was  noted  that  Preston  proved  that  a  rotating  ring  of  vortices  simulated  correctly  the 
energy  transfer  in  a  radial  flow  geometry,  the  details  were  not  examined.  Therefore,  it  is 
appropriate  now  to  look  further  into  the  application  of  vortex  loop  concepts  to  more  gen¬ 
eral  turbotnachines  having  centrifugal  effects. 

Before  proceeding,  the  term  “unsteady  potential  flow”  deserves  some  clarification. 
Examination  of  Eq.  (1)  shows  that  solutions  for  the  velocity  potential  depend  primarily 
on  geometric  boundary  conditions.  Thus  these  solutions  are  time-dependent  only  in  the 
sense  that  they  represent  instantaneous  velocity  fields  compatible  with  a  set  of  instanta¬ 
neous  boundary  conditions.  Therefore,  the  computed  velocity  fields  have  no  “memory”. 
The  pressure  field  is  different  since,  through  Eq.  (4),  it  depends  on  the  rate  of  change  of 
the  vector  velocity  and  will  affect  the  solution  through  any  boundary  condition  related  to 
pressure. 

For  generallized  turbotnachines,  the  standard  kinematic  botmdary  condition  applied 
to  Eq.  (1)  is  that  the  relative  velocity  normal  to  any  solid  surface,  stationary  or  moving,  is 
zero  (i.e.,  the  relative  velocity  is  always  tangent  to  the  surface).  This  condition  represents 
compensation  for  any  local  angle  of  attack  on  the  surface  by  the  undisturbed  flow.  For 
airfoil  shapes  with  a  sharp  trailing  edge,  the  boundary  condition  at  that  point,  that  substi¬ 
tutes  for  viscosity,  is  that  the  pressure  be  continuous  and  finite,  (e.g.,  a  stagnation  point.) 
In  matching  this  condition  with  vortex  singularities  it  is  necessary  to  locate  vortex  lines 
along  the  blade,  with  the  shed  vortex  system  convected  along  relative  flow  streamlines, 
such  that  the  flows  from  opposite  sides  merge  smoothly  at  the  trailing  edge. 
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Since  the  relative  velocity  is  so  critical  to  a  rotating  machine,  it  is  common  to  per¬ 
form  its  analysis  in  a  frame  of  reference  that  rotates  with  the  shaft.  Thus  the  relative  fluid 
velocity  in  a  frame  rotating  at  speed  Q  is  given  by; 

W  =  V-QxT  (29a) 

Because  the  inviscid  flow  is  assumed  to  remain  iirotational  in  the  absolute  frame,  i.e., 

^  s  0,  it  is  straightforward  to  show  that, 

Vxi^  =  -2Q  (29b) 

and,  in  general,  the  relative  velocity  will  be  rotational  and  not  subject  to  a  simple  poten¬ 
tial  flow  treatment.  Nevertheless,  W  can  often  be  treated  as  steady,  which  simplifies  the 
momentum  equation.  Making  the  strictly  kinematic  transformation  of  the  linear  momen¬ 
tum  equation,  (Eq.  (4),  into  a  rotating  frame  yields; 

O  «  ilK  =  zivp-V(H^V2)-(V  X  X  W-Q  X  (Q  X T)  +  2i^  X  Q  (30a) 
dt  Qo 

Substituting  (29b)  into  (30a)  shows  that  the  rotational  and  Coriolis  terms  cancel  to  give 

t\W  ^ 

o  » - =  V(p/o  -I-  W^/2)  -  Q  X  (Q  X  r)  =  ViRothalpy)  (30b) 

dt 

Equation  (30b)  is  therefore  a  rather  general  expression  for  computing  the  time-average 
local  pressure  in  the  rotating  frame. 

If  the  velocity  potential  in  the  absolute  frame  is  such  that  the  relative  flow  is  steady, 
then  it  can  be  seen  that  the  flow  field  rotates  with  the  blades,  and  (f>  (r,  0,  r)  becomes 
^  (r,0  -  Qr) .  Using  this  form  in  Eq.  (4)  gives; 


Rg.  8.  Illustration  of  “Coriolis'  force  in  the  rotating  frame. 
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and  for  steady  flow  in  the  rotating  frame,  Eq.  (30b)  beomes: 

o  =  -  V(p/e  +  -  {QrflD  ;  or  ViMj)  =  ViV^/2  +  iQr)^/2  -  W-/2)  (32) 

In  either  case  the  mean  change  in  absolute  frame  total  enthalpy  is  given  by  the  mean 
change  in  U$Vg  .  This  is  the  Euler  3-D  generallization  of  Eq.  (7).  If  the  flow  is  unsteady 
in  the  relative  frame,  then  the  relative  velocity  must  be  integrated  along  desired  paths  at 
successive  instants,  and  the  time  rate  of  change  computed  in  order  to  obtain  a  local  pres¬ 
sure  from  (30b). 

In  applying  the  vortex  singularity  model  to  a  purely  centrifugal  rotor,  it  should  be 
noticed  first  of  all  that  there  will  be  some  r~6  blade  profile  such  that  the  condition  of 
zero  relative  flow  normal  to  the  blade  surface  is  met  with  no  vortex  systems  required. 

This  “no-lift”  blade  shape  can  be  seen  to  have  a  greater  slope  as  the  radius  increases  due 
to  Eq.  (29a).  As  the  blade  shape  departs  from  the  no-lift  profile,  bound  and  shed  vortex 
lines  must  be  introduced  to  keep  the  flow  tangent  to  the  surface  and  meet  the  trailing  edge 
condition.  As  the  blade  becomes  primarily  radial,  the  bound  vortex  distribution  will  be 
strengthened  toward  the  outer  radii  because  of  the  higher  circumferential  velocities  to  be 
cancelled.  For  closely-spaced  radial  blades,  the  fluid  at  mid-radii  can  experience  flow  that 
is  mainly  radial  such  that  Wo  is  small.  In  this  region,  since  ^-0  and  the  rothalpy  is  con¬ 
stant  throughout,  the  following  approximations  apply; 


\v  = 


1  drWe 
r  dr 


1  dWr 
r  36 


-2Q;  -2eW,Ue 


(33) 


In  this  example,  significant  circumferential  variations  in  radial  velocity,  as  well  as  a 
Coriolis-like  distribution  in  static  pressure  occur  as  illustrated  in  Fig.  8.  Near  the  trailing 
edge  of  the  radial  blade,  relaxation  of  this  pressure  gradient  into  the  fluid  beyond  the  ro¬ 
tor  creates  a  flow  in  the  negative  direction.  This  reduces  machine  work  and  sets  up 
the  classic  “jet-wake”  behavior  described  by  Dean  and  Senoo  for  centrifugal  rotors.  Cer¬ 
tainly,  the  overall  behavior  is  described  by  the  system  of  rotating  vortices,  which  when 
summed  for  each  blade,  act  as  described  by  Preston.  But,  the  unique  pressure  distribution 
for  centrifugal  turbomachines  causes  a  non-uniform  exit  flow  field  that  will  be  shown 
later  to  be  acoustically  important  when  it  encounters  a  downstream  stator. 

ACOUSTICS  OF  DIPOLES 
PROPAGATION  OF  THE  VELOCITY  POTENTIAL 

For  various  kinds  of  mrbomachines,  the  blades,  rotating  or  fixed,  have  been  shown 
to  result  in  the  introduction  of  vortex  loops  or  equivalently  dipoles,  into  the  flow  field. 

For  rotating  blades,  the  time-varying  velocity  potentials  from  these  singularities  account 
for  the  power  transfer  to  the  fluid,  while  if  the  blades  are  stationary  the  energy  level  of  the 
fluid  is  constant  and  strictly  Bernoulli  changes  occur.  'Riming  to  the  acoustics,  it  is  now 
appropriate  to  ask  how  these  primarily  incompressible  potential  flow  effects  are  trans¬ 
formed  into  acoustic  signals  that  propagate  over  large  distances.  To  introduce  the  acoustic 
analysis,  it  will  be  assumed  that  the  fluid  density  can  in  fact  experience  very  small  varia¬ 
tions  with  pressure  through  an  isentropic  equation  of  state  such  that; 

dg  =  dp/c^  (34) 
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Writing  the  equations  of  conservation  of  mass  and  momentum  for  the  absolute  frame 
gives 

^+VqV  =  0  (35) 

dt 

^  =  -—V(p  +  gV^/2)  »  — —Vp  (36) 

ar  e  Qo 

Here,  following  Lamb,  the  kinetic  energy  part  of  the  total  pressure  can  often  be  neglected 
for  acoustic  purposes.  Using  (34)  in  (36)  and  integrating  in  time  gives  a  relationship  be¬ 
tween  velocity  and  density  perturbations. 

V  =  =  -  V  J piQodt  j 

Using  this  result  in  Eq.  (35),  with  the  condition  that  \V/c\,  the  local  Mach  number,  is 
small  gives 

(38, 

Equation  (38)  is  the  basic  wave  equation  having  solutions  made  up  of  functions  that  prop¬ 
agate  in  space  at  velocity  c.  Using  the  same  approach  it  is  possible  to  show  that  other 

variables  such  as  p,  V  or  6q  will  also  propagate  as  functioits  delayed  by  the  time 

t  =  j:/c. 

DEFnsrmoN  of  dipole-impulse  acoustics 

For  acoustic  analysis  it  is  often  convenient  to  locate  any  unsteady  activity  near  the 
origin  of  the  coordinate  system  and  to  determine  the  value  of  the  propagating  function  at 
some  remote  observation  point.  For  a  three-dimensional  spherical  geometry,  with  a  vol¬ 
ume  source  at  the  origin,  it  is  straight  forward  to  show  that  the  product  ^r  will  propagate 
radially  outward  from  the  origin  with  constant  amphtude.  The  velocity  potential  for  the 
dipole  is  constructed,  therefore,  by  combining  the  solutions  for  two  time-varying  volume 
sources  of  opposite  sign,  located  near  the  origin,  and  separated  by  a  small  distance,  2a, 
along  some  time-varying  direction  vector,  /tp  .  If  the  strength  of  the  dipole  is  defmed  as 

D{t)  =  Lim  [2aQ{t)]  (39) 

a-^o 

then  the  velocity  potential  at  any  point  in  the  field  of  Fig.  9  is  given  by  Eq.  (40)  which 
includes  behavior  both  near  and  far  from  the  migin. 


,  -  r/ D{t) 


Looking  at  the  0r  product  it  can  be  seen  that  the  grouping 
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will  propagate  to  large  distances  unattenuated.  As  is  conventional  in  acoustic  terminolo¬ 
gy,  the  square  brackets  indicate  that  the  function  inside  is  to  be  evaluated  at  the  retarded 
time. 


Fig.  9.  Geometry  for  dipole  radiation. 


Therefore  the  acoustic  velocity  potential,  and  all  the  other  far-field  properties  that 
can  be  derived  from  it,  depends  on  the  time  derivative  of  the  source  strength.  Using  the 
equivalencies  of  Eq.  (25)  it  can  be  seen  that  for  a  dipole  that  simulates  a  vortex  loop. 


Pvair,  t)  =  -Qo 


^va 

dt 


(42a) 


and  for  the  more  general  case  described  by  Lighthill  in  Eq.  (28)  for  a  body  moving  in  a 
fluid,  the  acoustic  pressure  is; 


(42b) 


Equations  (42)  are  the  ubiquitous  acoustic  results  for  the  application  of  a  force  to  a  small 
fluid  region.  For  a  turbomachine  with  multiple  blades,  the  acoustic  pressure  from  each 
source  is  summed  to  obtain  the  value  at  any  observation  point.  Although  the  procedure  of 
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applying  Eq.  (42)  appears  simple,  it  is  in  fact  difficult  because  of  the  critical  importance 
of  the  time  delay  operation  when  summing  signals  that  originate  from  different  parts  of 
the  flow.  Before  proceeding  with  techniques  for  the  evaluation  of  Eq.  (42)  for  turboma¬ 
chinery,  it  wiU  be  revealing  to  examine  first  the  simple  classic  case  of  dipole  acoustics 

VELOCITY  POTENTIAL  FOR  THE  OSCILLATING  SPHERE 

As  developed  previously,  the  dipole  is  related  to  the  vortex  loop,  and  to  the  impulse 
that  creates  such  a  flow.  In  addition,  ^e  dipole  singularity  can  represent  other  effects  such 
as  the  How  field  associated  with  the  motion  of  a  rigid  sphere  or  the  pressme  on  a  surface 
caused  by  a  time-varying  force  within  the  surface.  Considering  only  the  sphere,  there  are 
several  cases  of  acoustic  interest,  e.g.; 

•  Stationary  sphere  in  a  time-varying  flow 

•  Moving  sphere  in  a  stationary  fluid 

•  Sphere  oscillating  at  a  small  amplitude 

•  Sphere  rotating  in  a  plane  about  the  origin 

In  each  case  it  is  necessary  to  compute  the  remote  acoustic  velocity  potential  from 
some  local  boimdary  conditions  near  the  origin.  In  general,  this  problem  in  potential  flow 
is  handled  using  the  Green’s  integral  method,  including  the  time  delay  operation.  Thus,  if 
the  potential  is  known  on  some  surface  in  a  three-dimensional  flow  field,  the  potential  at 
some  remote  point  is; 


where  in  the  second  integral  only  the  time  delay  part  of  the  surface  derivative  is  usedThis 
equation  is  generally  difficult  to  evaluate,  except  in  simple  cases  where  the  surface  poten¬ 
tial  is  known.  One  of  those  cases  that  will  serve  as  an  example  is  that  of  the  rigid  sphere 
oscillating  in  a  fixed  direction  at  small  amplitude,  as  shown  in  Fig.  10.  For  this  case  a  di¬ 
pole  potential  can  be  defmed  easily  and  far-field  solutions  obtained.  To  do  this  it  is  first 
assumed  that  the  scalar  dipole  strength  can  be  described  as  a  periodic  function  that  is  in¬ 
dependent  of  direction. 

D(t)  =  Djijo))e^‘  (44) 

This  must  satisfy  the  boundary  condition  at  the  rigid  surface  that: 

Urir,fi,t)  =  —  \  r=R  =  Uz^e^' cos P  (45) 


Since  for  any  observation  point,  ■  ^o  =  cos  ^ ,  it  is  possible  to  combine  Eqs.  (40), 
(44)  and  (45)  to  obtain  an  expression  for  Do(Jq)). 


Do(j(0)  =  TjcR^U-,^ 


(l-ikR)^/2-jkR) 

{{l-kRf/2f+{kRf 


(46) 
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ng.  10.  Sphere  oscillating  at  small  amplitude  in  Z  direction. 


Here  k  -  o)/c  =  2n/k  and  it  is  generally  assumed  that  the  acoustic  wavelength  will  great¬ 
ly  exceed  the  sphere  circumference,  i.e.,  <  1.  Thus  for  this  case  the  complete 

velocity  potential  is: 


t) 


Do(Jo})  ^ja)Do(j(o) 


Ato-c 


cos^ 


(47) 


Putting  (47)  into  (43),  where  the  surface  of  integration  is  chosen  to  be  that  of  the 
sphere,  gives  for  the  far  field;  (here,  r  -  r^-Rno  ■  and  only  the  near  field  pan  of  (47) 
is  used) 


<pa(.r,t)  =  - 


cos^ 

AtVoC 


^  dU,  A  dU. 
dt  2  dt 


h-rfc 


(48) 


Thus  the  acoustic  velocity  potential  is  composed  of  two  terms,  one  depending  on  the  ac¬ 
celeration  of  the  displaced  volume  and  the  other  on  the  acceleration  of  half  the  displaced 
volume,  or  the  “added”  volume.  This  result  is  general  and  is  obtained  in  the  same  way  for 
any  propagating  scalar.  Although  (48)  could  have  come  from  (46)  and  (47),  the  use  of 
(43)  provides  the  acoustic  parallel  to  Eq.  (25),  and  the  analogy  between  (28)  and  (48). 
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POWER  INPUT  TO  THE  OSCILLATING  SPHERE 


Having  Eq.  (47),  it  is  of  interest  from  the  energy  viewpoint  to  examine  the  power 
flow.  Since  no  work  is  done,  other  than  acoustic  radiation,  this  wUl  give  some  insight  into 
the  source  of  the  acoustic  power  alone.  Looking  at  the  sphere  itself,  the  time-average 
power  input  is  determined  from  the  product  of  force  and  velocity  in  the  z  direction. 

/’w  =  ^  Re  [F,ijto)  •  U:*  Ob)}  (49) 

The  time-varying  force  in  the  z  direction  is  found  by  integrating  the  static  pressure  over 
the  surface  of  the  sphere. 


=r 


cosfi  2jt  R&infi  R  d  ^ 


From  Eq.  (5)  the  static  pressure  is; 


t)  =  -Qo^  -  Qoil^r  +  W 


The  component  surface  velocities  can  be  computed  from  (47)  as 


Ur  =  -^\  ^=U,cosfie^‘ 

dr  r=R  -o  ^ 


1  d<t>  .  <t>{.R,t) 


tan)S  (52) 


It  can  be  seen  that  the  angular  dependencies  of  the  kinetic  energy  terms  are  such  that 
when  used  in  (50)  the  result  will  be  zero.  Thus  the  z  direction  force  becomes  only  a  func¬ 
tion  of  the  first  term  in  (51). 

Fzit)  =  J  QoDo(j(o)j(0  ^  2n  R^  sinfidfi  (53) 


The  complex  value  of  the  integral,  using  (46)  for  Do(jo>) ,  and  assuming  kR  is  small,  is 

Fzijo))  =  Qo^U;Q}(j  +  [kRf/2)  (54) 

Equation  (54)  shows  that  there  is  an  inertial  force  on  the  sphere  associated  with  the  added 
volume,  A/2  with  a  small  fraction  of  that  force  being  in-phase  with  velocity.  Therefore 
from  (49)  the  mean  power  input  needed  to  oscillate  the  sphere  is: 


2 


POWER  PROPAGATING  FROM  SPHERE 

To  obtain  this  result  from  a  different  view,  one  can  also  compute  the  acoustic  power 
that  radiates  to  infinity  from  the  origin.  Thus  in  the  far  field; 
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r  ,x___  ^  _  -gJatPoUQ})  j(0  ^ 

r,  f)  —  Qq  COS^ 


"  "  dt  4jrr  c 

,  ..  ^  -joiDoijoi)  jo) 

Ua(.r,t)  =  -Qo—  = - 5 - COS^  (56b 

dr  4nr  cr 

If  can  be  seen  that  p.  =  g^u, ,  as  in  plane  wave  propagation,  and  for  kR«l  both  coeffi¬ 
cients  are  real  and  in-phase.  Their  product  gives  the  time  average  acoustic  intensity. 

_  gc0  Dli^  2  a 

laiTo^p)  =  - 75-cos^^  (57) 

KAjtToY 

Integrating  this  over  an  arbitrarily  large  spherical  surface,  with  Do  =  TsiR^Uzo  gives; 

1  /  AtcR^  \  ( U,\  (krY 


Pout 


1  /  AtcR^  \  ( U,\ 


Which  is  exactly  equal  to  the  input  power,  Eq.  (55).  Looking  at  (56a)  more  closely. 


r»,/!)  =  3(M) 


U,Q)  Q) 

^TCToC 


which  suggests  from  Eq.  (42b)  that  the  force  involved  in  creating  the  acoustic  signal  is 
proportional  to  the  inertia]  acceleration  of  the  added  mass  and  the  displaced  mass  terms. 
Thus  the  total  apparent  force  experienced  by  the  fluid  is  found  by  considering  the  inside 
of  the  sphere  as  being  filled  with  the  surrounding  fluid.  This  is  consistent  with  Lighthill's 
expression  for  the  total  force,  Eq.  (28),  and  the  general  solution,  Eq.  (48).  However,  now 
it  can  be  seen,  by  comparing  (28)  and  (48),  that  the  time  derivative  of  the  volume  integral 
of  the  vorticity  moment  must  be  equal  to  the  acceleration  of  the  added  mass.  This  is  in 
fact  the  case  because  the  displaced  mass  effect  is  associated  with  the  radial  velocity  of  the 
sphere  while  the  added  mass  is  associated  with  the  tangential  velocity  on  the  sphere.  The 
latter  effect  is  equivalent  to  a  set  of  vortex  loops  in  the  surface  of  the  sphere  that  account 
for  the  slip  on  this  surface.  The  combined  effect  of  this  set  of  ring  vortices  is  one  single 
vortex  loop  with  its  axis  in  the  z  direction.  In  the  far-field  this  vortex  loop  appears  as  a 
single  dipole  representing  the  effect  of  the  added  mass.  Through  Eq.  (54),  tMs  added 
mass  dipole  includes  a  small  component  in-phase  with  the  velocity,  and  therefore  ac¬ 
counts  for  all  the  acoustic  power  input. 

To  summarize  the  acoustic  behavior  of  the  oscillating  sphere,  the  total  force  acting 
on  the  fluid  as  a  result  of  the  motion,  and  which  is  the  apparent  source  of  radiation  into 
the  far  field,  is  due  to  the  inertia  of  both  the  "displaced”  and  “added”  masses.  Interesting¬ 
ly,  the  displaced  mass  term  cannot  contribute  to  the  acoustic  power  since  its  force  is 
always  in  quadrature  with  the  velocity.  However,  it  so  happens  that  the  added  mass  term 
contains  a  small  component  in  phase  with  velocity,  which  accounts  fully  for  the  acoustic 
power.  In  addition,  the  added  mass  term  is  related  to  a  disnibution  of  vortex  rings  on  the 
surface  of  the  sphere,  and  thus  to  a  single  dipole,  and  fmally  to  the  impulse  needed  to  set 
up  an  external  tangential  flow  around  the  sphere.  Thus,  although  it  can  be  seen  from  Eq. 
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(40)  that  the  tangential  velocity  itself  does  not  propagate  to  the  far-field,  the  impulse  re¬ 
sponsible  for  it  does.  More  to  the  point,  this  particular  vortex  ring  result  is  related  to  the 
vortex  loops  introduced  by  the  work  input  action  of  turbomachine  blades.  Thus,  In  the 
case  of  the  work  vortices,  only  Eq.  (42a)  is  necessary  if  the  blade  volume  is  negligible. 
Finally,  when  extrapolating  from  the  near  field  to  the  far  field  it  is  important  to  realize 
that  there  are  two  components  to  Eqs.  (28)  or  (43).  One  author,  Blake,  used  only  the  ra¬ 
dial  surface  velocity  to  compute  far  field  pressure,  and  obtained  two-thirds  of  the  answer. 
The  use  of  a  scalar  function  instead  of  a  vector  to  extrapolate  to  the  far-field  will  prevent 
this  type  of  error. 

ACOUSTIC  INTERFERENCE  FOR  TURBOMACHINES 

Turbomachines  represent  a  unique  case  in  the  computation  of  the  far-field  sound 
pressure  from  Eq.  (42)  because  the  dipole  sources  are,  by  definition,  located  in  a  circular 
pattern  about  the  rotational  axis.  If  the  sources  are  the  steady  work  forces  that  the  rotor 
blades  apply  to  the  fluid,  then  they  simply  rotate  in  space  at  constant  amplimde.  On  the 
other  hand  if  the  acoustic  sources  are  the  dynamic  forces  on  stator  blades,  perhaps  due  to 
the  wake  of  an  upstream  rotor,  then  the  sources  are  stationary  in  space,  while  the  exitation 
pattern  rotates.  In  all  cases  the  common  requirement  for  solving  Eq.  (42)  is  to  determine 
how  time  delays  can  be  accounted  for  with  a  circular  array  of  sources.  To  illustrate  the 
approach  to  this  problem  for  low  Mach  number  applications,  two  basic  acoustic  propaga¬ 
tion  geometries  will  be  examined  for  steady  shaft  rotation  about  a  point  fixed  in  space; 

•  Radiation  into  a  free-fleld 

•  Radiation  in  a  cylindrical  duct 

For  each  case  the  solution  will  be  examined  for  a  rotor  alone,  for  a  stator  behind  a 
rotor,  and  for  a  rotor  operating  in  a  spatially  non-unifoim  inflow.  In  cases  where  the  flow 
exitation  is  periodic,  the  force  response  will  be  assumed  to  be  linear  in  the  perturbed  flow. 
For  an  isolated  blade  or  a  cascade,  the  appropriate  Sears  function  will  give  the  unsteady 
force  for  the  airfoil,  while  the  “Coriolis”  force  must  be  used  for  rotors  with  radial  flow 
conditions.  It  will  also  be  assumed  that  the  effects  of  fmite  blade  volume  can  be  described 
by  the  potential  source  functions  described  previously  for  the  solid  sphere  as  in  (Quandt 

FREE-FIELD 

Gutin  Sound  From  a  Rotor 

This  case,  representing  a  stationary  open  propeller,  was  first  analyzed  correctly  by 
Gutin,  and  the  res^t  is  termed  “Gutin  Sound”.  Because  the  solution  is  basic  to  all  later 
analyses,  the  approach  will  be  described  in  some  detail.  For  steady  inflow  to  a  single  pro¬ 
peller  blade,  the  work  force  is  seen  to  be  constant  in  magnimde,  but  varying  in  direction 
and  point  of  origin  as  the  blade  rotates.  Since  this  problem  is  symmetric  about  the  z  axis 
of  rotation,  it  is  convenient  to  locate  the  observer  in  the  x-z  plane,  as  shown  in  Fig.  1 1 . 
Thus  the  blade  force  vector  becomes 

F  =  I  F  I  ns  =  F{-T siny  sin5,+y  siny cos0,-^'cos  y)  (60) 

where  y  is  the  angle  between  the  blade  zero  lift  line  and  the  x  axis.  The  unit  vector  to  the 
observer’s  position  is; 


DTRC-PAS-91/55 


21 


Fig.  11.  Acoustic  geometry  for  a  single  blade  rotating  in  the  x-y  plane. 

no  =  ismficosip  +  kcosfi  (61) 

and  the  function  that  the  observer  in  the  hemisphere  {if>  =  0)  will  see  is: 

[no-^],^/^  =  -^[-/"sinysm^sine,  +  Fcosycos/S]^^^^  (62) 

In  Eq.  (62)  the  second  term  represents  an  axial  force  that  is  steady  in  time  and  direc¬ 
tion,  but  varying  in  origin  as  the  blade  rotates.  It’s  time  derivative  gives  rise  to  a  term  that 
is  linear  in  Mach  number,  and  thus  it  will  be  neglected  at  this  point  The  first  term  is  more 

clearly  time  dependent  since  d,-Qt.  It  represents  the  torque  force  on  the  blade,  and  as 
such  can  be  seen  to  have  force  components  that  oscillate  in  the  x  and  y  axes.  Gt  may  be 
noted  that  the  observer  in  the  x  -  z  plane  will  see  no  variation  in  signal  from  the  compo¬ 
nent  along  the  y  direction). 

With  a  force  that  repeats  periodically  at  frequency  Q,  it  can  be  assumed  that  the 
observer  will  also  experience  an  acoustic  pressme  from  Eq.  (42)  that  is  periodic  at  Q, 
such  that; 

oe 

Pa{r,P,t)=  X  A„{r,P)  (63) 

ns- 00 
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and  it  is  desired  to  determine  the  A„  terms.  Now,  any  signal  reaching  the  observer  at  time 
t  had  to  originate  at  some  earlier  time.  By  the  geometry  of  the  system,  the  distance  from  a 
blade-fixed  source  rotating  at  some  radius  about  the  origin,  to  the  chosen  far  field  observ- 


(  ^  ■ 
=  ro|  1 - SI 

V 


sinfi  cosdg 


Therefore  the  angular  position  of  the  blade  at  an  earUer  time,  determined  by  the  acoustic 
delay,  at  which  a  signal  originated  that  reaches  the  observer  at  t  is: 


[e,]  =  Q  ^r- 


Solving  Eq.  (65)  for  «  Q  r  and  substituting  in  Eq.  (63)  gives; 


=  n-  [0,]+ - 

c 


|cos[0J-0, 


where  for  simplicity  the  following  definition  has  been  used: 

CIR 

I  = - sin  =  Mj-sin  (67) 

c 

Combining  Eq.  (42)  with  the  general  expression  for  the  coefficients  in  Eq.  (63)  gives 

L„ei’>^Qr,/c~e.o)  L  f  (6g) 

^jCToC  2  dt 

In  order  to  evaluate  this  Fourier  integral  it  is  most  advantageous  to  use  the  following 
identity  uniquely  appropriate  to  circular  sources,  as  introduced  by  Gutin. 

00 

(69) 


Using  (69)  in  (68)  and  interchanging  the  order  of  summation  and  integration  allows  the 

Fourier  integral  to  be  written  as: 

1  _  .-ae.+ne.-na.) 

T-  - 7^ -  dO,  (70) 

2Jil„  2  j 


This  integral  has  a  value  of  1/2  only  when  m  =  n  ±  1,  otherwise  it  is  zero.  Therefore  in  the 
evaluation  of  An ,  there  appear  two  Bessel  functions,  which  can  be  combined  to  give  a 
single  Bessel  function  of  order  n.  The  resulting  far  field  acoustic  pressure  for  one  blade 
becomes; 

=  T  (71) 
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If  now  the  turbomachine  is  composed  of  B  equally  spaced,  and  equally  loaded  blades,  the 
acoustic  pressure  is  the  sum  over  the  individual  Bso  positions.  Since  this  sum  is  B  for 
n  =  lB  ,  and  zero  otherwise,  the  result  is; 

PaBir,^,  t)  =  X  (72) 

which  is  the  classic  Gutin  solution  to  Eq.  (42)  for  low  Mj. 

Some  illustrative  features  of  the  acoustics  of  a  single  rotor  in  a  free  acoustic  field 
can  be  obtained  by  examining  Eq.  (72).  First  of  all,  the  leading  term  is  the  total  force  on 
the  rotor  blades,  radiating  in  such  a  manner  as  to  have  the  acoustic  pressure  decrease  in¬ 
versely  withr^  .  Secondly,  the  amplitudes  of  the  steady  term,  /  =  0  ,  the  fimdamental,  and 
the  higher  frequency  terms  are  all  indexed  to  the  Bl  product,  with  a  spatial  interference 
pattern  given  by  the  Bessel  function  of  order  Bl . 

To  interpret  this  result  it  is  useful  to  know  how  Bessel  functions  behave  versus  their 
argument,  as  shown  in  Fig.  12  forJ^,  Ji  mdJ^ .  The  following  asymptotic  relationships 
are  useful: 


»  1.0 

for^  =  «  1 

(73a) 

1  JaiiBli)  1  « 

Bl 

(73b) 

f  2 

1  JBKBm  1  “ 

(73c) 

Fig.  12.  Representative  Bessel  functions  vs.  argument. 
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Now,  when  /  =  o,  the  result  is  always  zero  since  this  index  value  multiplies  each  harmon¬ 
ic.  As  the  harmonic  index  increases,  with  ^  still  small,  the  amplitudes  fall  rapidly  with 
increasing  /.  In  all  cases  the  acoustic  pressure  on  the  axis,  ^  =  o  will  be  zero.  Since  the 
Bessel  functions  oscillate  about  zero  as  the  argument  increases,  the  acoustic  pressure  ex¬ 
hibits  a  symmetric  lobed  pattern  with  increasing  angle  ^  .  In  general,  as  the  harmonic 

number  increases,  the  amplitudes  of  the  pressure  pattern  will  decrease  quickly.  The  impli¬ 
cation  of  the  basic  Gutin  solution  is  that  larger  blade  numbers  should  reduce  the  tonal 
acoustic  pressure  for  the  low  tip  Mach  numbers  of  interest. 

Stator  In  the  Wake  of  a  Rotor 

The  second  arrangement  to  examine  is  that  for  a  stator  vane  row  situated  down¬ 
stream  of  a  rotor.  A  convenient  technique  for  this  geometry  was  reported  by  Duncan  and 
Dawson,  who  were  concerned  primarily  with  the  total  acoustic  power  radiated  from  un¬ 
equally  spaced  blades,  instead  of  a  detailed  description  of  the  radiated  acoustic  pattern. 
Duncan  and  Dawson  assumed  an  unsteady  forced  i.e.,  dipole,  radiation  from  each  stator, 
and  that  the  force  on  any  stator  blade,  q,  was  attributable  to  the  wake  from  rotor  blade,  b. 
Figure  13  shows  how  the  motion  of  each  blade  produces  a  variation  in  force  on  each  sta¬ 
tor  as: 

=  i  (74) 

l=~a> 

Here  the  Fourier  coefficient  is 

Bi  =  F(0  =  f  ddB  (75) 

Where ffiBa)  is  the  single  event  blade-vane  interaction  centered  a\6  =  o  and r  =  o  ,  and 

the  event  is  repeated  for  every  ^6  =  2ji  .  If  now  there  are  B  equal  and  equally  spaced 
blades,  the  summed  force  on  any  stator  vane  is; 
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(76) 


00  B-\ 

Fqi{eB)=  X 

/=-oo  b=o 

Since  df,  ~  htbjB  and  Ob  =  Q  t  ,  the  sum  is  non-zero  only  for  /  =  nB  such  that; 

00 

F^i{t)  =  B  X  ill) 

n=-oo 

gives  the  force,  and  therefore  the  dipole  strength  on  each  stator  vane.  From  (42)  it  is 
known  that  this  force  is  radiated  widi  the  acoustic  delay  time; 

,  R  simB 

tq  =  Tote - ^-^cos(V'  -  (78) 

Therefore  the  far-field  acoustic  pressure  is: 

Pair,  t)  =  X  (79) 

where  ^  is  defined  in  Eq.  (67).  The  direction  vector  product  is; 

'nonq  =  cosy  cos/S  +  siny  sin^  sm(V'  -0^)  (80) 

which  again  represents  axial  and  tangential  force  components.  Introducing  the  Gutin  Bes¬ 
sel  function  identity  from  Eq.  (69)  and  summing  over  all  V  equal  vanes  gives; 

Pa(r,t)  =  -j^  X  X  X 

ns-co  ms-se  q=o  (81) 

Comparing  the  overall  features  of  (81)  and  (72),  it  can  be  seen  again  that  the  acoustic 
pressure  varies  inversely  with  rg  ,  and  that  the  Bessel  function  behavior  is  symmetric 
about  the  plane  of  rotation. 

Inspecting  (81)  in  more  detail,  it  can  be  seen  that  the  sum  over  all  vanes  for  the  axial 
term,  /3  =  o,  of  (80)  will  be  zero,  except  for; 

m  +  nB=±kV  (82) 

In  this  case,  when  the  blade  and  vane  numbers  are  equal,  or  they  have  some  common 
multiple,  the  index  m  =  o  occurs,  and  the  Jo  term  can  appear.  This  corresponds  to  a  con¬ 
dition  where  all  the  blade  wakes  and  vanes  overlap  axially  at  the  same  time  such  that 
there  is  unified  dipole  radiation  in  the  axial  direction.  Compared  with  the  rotor  alone,  this 
radiation  can  be  quite  strong,  even  at  low  Mach  numbers,  because  of  (73a).  However  it 
has  been  known  for  some  time  that  this  tonal  soiuce  can  be  eliminated  by  making  V  >  2B 
such  that  these  Jo  terms  do  not  develop.  Of  course,  the  strength  of  each  harmonic  will 

vary  with  sharpness  of  the  wake  pulse  coefficient,  B„b  ,  which  is  in  turn  related  to  a 
Sears  function  response. 
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Turning  to  the  torque  term,  Eq.  (80)  the  indexed  position  of  each  source  relative  to 
an  observer  in  the  x-z  plane,  i.e.,  when  \p  -0  introduces  conditions  such  that  the  sum 
over  all  vanes  is  non-zero  when: 

m  =  nB-^kV+\  (83) 

Again  there  is  the  possibility  for  m  =  o,  and  Jo  to  appear  when  the  blade  and  vane 
number  combinations  are  separated  by  an  integer  of  one.  The  implication  here  is  that  the 
radiation  in  the  radial  direction  at  frequency  nSQ  behaves  as  a  rotating  dipole.  Again,  the 
harmonic  strength  depends  onfi„  ^  ,  and  the  Bessel  function  order  can  be  kept  large 
whenV  >  2B  . 

Rotor  in  a  Non-Untform  Inflow 

The  last  free-field  radiation  case  to  examine  is  that  for  a  rotor  operating  in  a  steady, 
but  non-uniform,  inflow  as  shown  in  Fig.  14.  The  inflow  distortion  can  be  described  by 
circiunferential  harmonics  of  the  form  of  (84)  if  it  is  deterministic: 

dvjs  =  Vh.  cos  J  (84) 

When  this  distribution  modulates  the  steady  blade  force  used  in  the  Fourier  integral  of 
(68),  the  sindf  term  from  the  steady  force  acquires  a  perturbation  component 

sin  (h'-h  1)05--^  sin(w- 1)0,^  (85) 

Using  this  in  Eq.  (70),  it  can  be  seen  that  the  only  Bessel  function  orders  having  non-zero 
coefficients  are: 


Rg.  14.  Rotor  blade  operating  in  a  non-uniform  deterministic  inflow. 
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m  =  w  +  (M’+l)  (86) 

m  =  n  +  iw—  1) 

By  combining  the  resulting  Bessel  functions,  and  summing  the  acoustic  s’'nals  from  B 
equally  spaced  blades,  the  far-field  pressure  is; 

=  i  (87) 

2 

From  (87)  it  can  be  seen  that  if  the  inflow  disturbance  number  is  some  multiple  of  B,  then 
Jo  terms  appear  and  strong  on-axis  dipole  radiation  occurs,  even  at  low  Mach  number. 
However,  even  if  B  and  w  have  no  common  multiples,  but  are  such  that  low  order  Bessel 
terms  appear,  i.e.,  Ji,  or  Jo,  then  some  signiBcant  acoustic  radiation  can  still  occur  at  blade 
passing  frequencies,  being  strongest  in  the  radial  direction. 


CYLINDRICAL  DUCT 

When  the  turbomachine  rotor  is  located  within  a  duct  system  it  is  necessary  to  ac¬ 
count  for  the  boundaries  in  computing  solutions  to  Eq.  (42).  For  internal  flow  the 
simplest  geometry  is  the  rotor  in  an  infrnitely  long  cylindrical  duct  with  rigid  walls.  This 
geometry  is  appropriate  for  axial  flow  machines,  but  for  centrifugal  devices  a  toroidal 
volute  is  often  present  between  inlet  and  exit  ducting.  That  geometry  produces  its  own 
unique  acoustic  cavity  behavior,  which  must  be  examined  separately  using,  for  example, 
the  finite  element  approach  of  Everstine.  In  the  present  study  the  focus  will  be  limited  to 
axial  flow  machines  in  a  cylindrical  duct. 


Tyler  and  Sofrin  Rotor-Stator  Tones 


The  classic  solution  to  the  acoustic  radiation  from  an  axial  flow  rotor,  or  rotor-stator 
combination,  in  a  duct  is  that  of  IV^cr  and  Sofrin.  Their  approach  considers  first  the  gen¬ 
eral  question  of  acoustic  propagation  of  unsteady  pressure  patterns  in  a  cylindrical  duct, 
rather  than  specific  blade  dipoles.  Thus  it  is  found  that  only  certain  circumferential  mode 
shapes,  associated  with  specific  radial  Bessel  function  patterns,  will  propagate  over  sig¬ 
nificant  axial  distances.  In  particular,  only  those  compatible  radial  and  circumferential 
pressure  distributions  that  appear  to  rotate  at  greater  than  sonic  sp>eed  are  found  to  propa¬ 
gate  axially.  Thus  a  rotor-locked  pressure  field,  i.e.,  steady  in  the  rotating  frame,  will 
propagate  only  if  the  rotor  speed  is  supersonic.  On  the  other  hand,  the  stator-locked  pres¬ 
sure  field,  which  behaves  as  a  standing  wave  system  excited  at  blade  passing  frequency, 
can  be  thought  of  as  two  oppositely  rotating  pressure  fields  since; 


cos(5Qr)  cos 


TjtRO  1 


( 

cos  I 


IstRB  \  /  ^  TjtRO 

BQt  +  — - I  +cos|  BQt - - — 
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(88a) 


Rewriting  the  second  term  gives,  for  example; 
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2jt  j 
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(88c) 


so  that  the  apparent  circumferential  Mach  number  for  a  pressure  wave  lengthAu’  is: 

BQ  Aw  Aw 

Mw  = - = - —  =  c — 

c  c  2n  kag 

Therefore  when  the  circumferential  wavelength  of  any  Fourier  pressure  component  ex¬ 
ceeds  the  acoustic  wavelength  for  a  frequency  corresponding  to  the  repeating  blade 
exitation  rate,  then  the  wave  speed  will  appear  to  be  supersonic.  The  result  is  that  smaller 
spatial  wavelength  disturbances  will  not  propagate. 

With  the  above  criteria  in  mind,  lyier  and  Sofrin  assume  that  the  rotor/stator  system 
of  Fig.  is  will  create  the  following  annulus  pressure  modes  for  the  interaction  of  the 
blade  row  with  a  single  vane  located  at  the  angular  origin; 

OD 

Pq{e,t)=  X  (89) 

n=-oo 

where  To  =  TjiIBQ  .  Further,  it  is  assumed  that  each  Fourier  coefficient  in  Eq 

be  described  by  a  second  combination  of  circumferential  harmonics; 

00 

A„(0)=  X 

l?J=-00 

If  this  arbitrary  pressure  pattern  is  now  summed  for  multiple  identical  stators  with  angular 
separation  of  A0  =  lijt/V  ,  and  time  delay  r  =  A6/Q  from  the  start  of  the  sequential 
events,  the  annulus  pressure  is: 


.  (89) can 

(90) 


Rg.  15.  Multibladed  rotor  acting  on  a  stator  vane  in  an  infinite  duct. 


DTRC-PAS-91/55 


29 


(91) 


^  ^  ^  g^m(e-9A0)+nSQ(/+9Afl/Q )) 

/is— »  m^oo 

Combining  terms  into  appropriate  summation  groups  gives 

P0.t)=  i  i 

n=-oo  ma-oe 


(92) 


From  Eq.  (92)  it  is  clear  that  the  sum  over  all  the  vanes  of  the  last  term  will  be  zero  im- 
less; 

m  =  nB  +  kV  (93) 

for  which  it  has  the  value  V.  Additionally,  from  the  mode  summation  exponent  it  can  be 
seen  that  the  apparent  speed  of  rotation  of  the  mth  mode  shape  is; 

nB 

V„  =  —  OR  (94) 

m 


where,  from  the  initial  acoustic  analysis,  that  mode  will  propagate  if  V„  >  c. 

Even  though  the  lyier  &  Sofrin  results  are  very  simple,  quite  a  bit  can  be  learned 
from  their  examination.  For  example,  Eq.  (94)  indicates  t^t  all  modes  satisfying  the  m  = 
o  condition  will  propagate,  effectively  as  plane  waves,  since  kw  is  innnite.  This  is  similar 
to  the  appearance  of  Jo  radiation  by  Eq.  (82)  for  the  free-field  problem.  This  mode  is  as¬ 
sociated  with  thrust  and  torque  fluctuations  on  the  stator  from  all  vanes  acting  together. 
Inspection  of  Eq.  (93)  shows  that  m  =  o  when  B  and  V  are  equal,  or  have  some  common 
multiple.  In  case  neither  of  the  above  criteria  is  met,  then  the  first  value  of  blade  harmon¬ 
ic  n  that  gives  m  ==  o  will  be  V,  and  the  BV  harmonic  will  radiate,  although  the 
amplitude  ,  may  be  small.  Continuing,  if  there  is  one  blade  that  differs  from  the  oth¬ 
ers,  then  B  =  1  is  appropriate  and  m  =  o  occurs  for  n  =  V  such  that  vane  frequency 
propagates.  On  the  other  hand  if  there  is  one  odd  vane,  then  V  =  1  and  m  =  o  occurs  for 
k  =  B  and  «  =  1  such  that  blade  frequency  propagates.  For  the  rotor-only  case,  V  =  o 
and  Eq.  (94)  indicates  that  with  m  =  rB  ,  the  rotor  tip  speed  must  be  supersonic  for  any 
radiation  to  propagate  in  the  cylindrical  duct.  Finally,  even  if  m  =  1,  2, 3, ...  propagation 
can  still  occur  if  V„  from  Eq.  (94)  is  greater  than  c.  In  such  cases  the  wave  fronts  are 
angled  to  the  duct  axis,  therefore  giving  rise  to  the  term  “spinning  modes." 

As  a  general  conclusion,  lyier  and  Sofrin  predict  that,  for  a  subsonic  rotor-stator 
combination  in  a  duct,  tonal  acoustic  radiation  can  be  prevented  to  a  large  degree  when  V 
>  2B.  This  r\ile  is  similar  to  the  free-field  result  of  Eq.  (82),  which  shows  that  the  V  >  2B 
criteria  will  prevent  Jo  terms,  and  their  strong  radiation.  It  also  provides  a  constraint  in  the 
free-field  that  only  higher  order  Bessel  terms,  with  their  lower  ampUtudes,  wiU  appear  for 
subsonic  turbomachines.  According  to  Cumpsty  this  result  has  been  known  and  used  for 
many  years  in  the  aircraft  literature. 

Goldstein’s  Classic  Analysis 

The  approach  of  lyier  and  Sofrin  is  one  way  to  represent  the  effect  of  multiple  di¬ 
poles  acting  as  a  phased  planar  array  in  a  rigid  cylindrical  duct.  The  success  of  their 
analysis  suggests  that  the  exact  nature  of  the  pressure  field,  either  near  or  far,  is  not  as 
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important  as  is  the  spatial  pattern  and  the  time  delays.  In  fact,  Tyler  and  Sofrin  are  vague 
about  just  where  the  pressure  pattern  is  defined  to  start,  do  not  directly  require  the  Gutin 
identity  of  Eq.  (69),  and  therefore  do  not  produce  Bessel  functions  as  such,  with  physical¬ 
ly  derived  dipole  force  coefficients.  Of  course  the  Bessel  function  for  the  duct  problem 
will  be  different  than  for  the  free  field,  the  pressure  wiU  not  vary  inversely  with  distance, 
and  the  distinction  between  axial  and  radial  directivity  will  not  exist. 

These  imcertainties  and  similarities  have  been  resolved  by  Goldstein  who  develops 
exact  analytic  solutions  for  the  rotor-statCH’  problem  in  finite  and  infinite  length  ducts.  Us¬ 
ing  traditional  separation  of  variables  techniques,  Goldstein  starts  fiom  local  blade  forces 
in  terms  of  torque  and  thrust,  and  expresses  the  radial  pressure  patterns  in  terms  of  Bessel 
functions.  Applying  the  boundary  condition  of  zero  radial  velocity  at  the  rigid  duct  wall 
requires  that  the  pressure  gradient  also  be  zero,  and  therefore  that  the^rsr  derivative  of 
each  Bessel  function  be  zero  at  the  wail.  Since  Bessel  functions  oscillate  radially,  there  is 
an  infmite  series  of  zeros  for  each  order.  These  radial  root  numbers  enter  into  the  axial 
propagation  function  such  that  there  wiU  be  rapid  decay  of  the  mode  when: 


(1-M2) 


(95) 


Here,  is  the  s  order  zero  of  the  first  derivative  of  Jm-  Goldstein  points  out 


that  tends  to  infmity  with  m  or  s  such  that  only  the  initial  roots  need  be  examined. 


Since  the  first  zero  turns  out  to  always  be  larger  than  m/R,  propagaticm  will  not  occur 
when 


Examining  the  Bessel  function  orders  m,  that  describe  the  force  distribution,  Gold¬ 
stein  shows  that  for  a  deterministic  inflow  distortion  of  circumferential  order  w  ,  the 
following  are  permissible; 

(97a) 


and  for  a  rotor-stator  system  of  V  vanes  the  Bessel  function  orders  are  again: 

m  =  nB  +  kV 

In  all  cases  there  will  be  propagation  along  the  duct  axis  when: 

a 


(-"'J 


I-A/2 


(97b) 


(98) 


Therefore  for  the  rotor  by  its^  in  a  uniform  flow,  i.e.,  V  =  o,  there  will  be  propagation  in 
the  duct  when 

+  (99) 


For  wake  or  vane  generated  force  oscillations  there  will  always  be  plane  wave  propaga¬ 
tion  and  thrust  oscillations  when  the  Jo  terms  occur.  When  higher  order  Bessel  function 
shapes  satisfy  Eq.  (98)  the  propagation  will  occur  in  the  form  of  “spinning  modes”  as  de¬ 
scribed  first  by  Tyler  and  Sofrin  and  illustrated  in  Fig.  16. 
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Fig.  16.  Propagation  of  plane  and  spinning  acoustic  waves  from 
rotor-excited  sources  in  a  duct. 

While  Goldstein  reaches  much  the  same  conclusions  as  lyier  and  Sofrin,  his  analy¬ 
sis  is  more  definitive  in  that  the  forces  responsible  for  the  acoustic  waves  are  identified. 
Thus  the  coefficients  of  each  Bessel  function  are  exactly  the  linear  coefficients  of  the 
blade  response  to  a  periodic  change  in  upstream  velocity.  Hence  aU  coefficients  depend 
upon  a  Sears  function  solution  for  the  blade  or  vane  row,  such  that  the  ampUtudes  will 
decrease  rapidly  with  reduced  frequency  greater  than  one. 

In  summary,  the  analysis  of  Goldstein  for  turbomachine  acoustics  in  a  duct  agrees 
with  lyier  and  Sofrin,  but  provides  the  necessary  details  on  the  dipole  nature  of  the 
acoustic  source.  For  the  rotor-stator  combination,  the  conclusion  remains  that  tonal  radi¬ 
ation  should  not  occur  for  V  >  2S  .  Machines  that  exhibit  tones,  even  though  this 
criterion  has  been  met,  are  believed  to  be  influenced  by  steady  inflow  distortions  or  by 
large  scale  turbulence.  As  alluded  to  earher,  die  major  tonal  difference  between  a  rotor- 
stator  combination  and  a  stator-rotor  combination  is  that  the  harmonic  thrust  coefficients 
are  indexed  in  nB  for  the  former  and  kV  for  the  latter,  but  in  all  cases  the  radiation  must 
occur  with  period  determined  by  the  moving  member.  Finally,  for  all  interference  analy¬ 
ses  in  systems  with  circular  symmetry  the  presence  of  Bessel  functions  describes  the 
interaction  of  discrete  and  perfect  geometries.  In  real  physical  systems  many  imperfect 
conditions  occur  which  tend  to  blur  these  results.  When  the  turbomachines  have  been  ^- 
signed  to  eliminate  first  and  second  order  acoustic  signals,  then  other  types  of 
propagation  analysis  may  be  more  appropriate  than  the  perfectly  harmonic  ones  described 
above. 
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VORTEX  CONCEITS  IN  REAL  TURBOMACHINES 


VORTICES  AND  PERFORMANCE 

The  previous  analysis  showed  how  the  flow  structure  of  unsteady  vortex  loops  is 
associated  with  the  work  output  and  the  tonal  acoustics  of  turbomachines  for  some  simple 
cases.  While  this  provided  an  insight  into  the  large  scale  flow  mechanisms  involved,  there 
remains  the  question  of  how  these  concepts  can  help  to  defme  the  details  of  real  machines 
where  the  fluid  has  rotationality,  viscosity  and  turbulence.  From  a  larger  perspective,  this 
complete  problem  is  usually  investigated  by  quantitative  flow  visualization  or  by  compu¬ 
tational  fluid  dynamics,  CFD.  With  the  advent  of  modem  computers,  the  latter  technique 
allows  the  simultaneous  solution  of  the  iull  mass,  momentum  and  energy  conservation 
equations  in  detail  by  specifying  only  the  behavior  of  the  fixed  and  moving  boundaries. 
The  continuing  work  of  Lakshminarayana  and  Thompson,  and  Rai,  and  co-workers  are 
representative  of  CFD  developments  directly  applicable  to  the  turbomachines  of  interest 
herein.  Although  these  diagnostic  techniques,  experimental  or  computational,  can  give 
more  realistic  results  than  a  potential  flow  analysis,  some  of  the  phenomena  unique  to 
turbomachine  performance  and  acoustics  can  still  be  appreciated  most  readily  in  terms  of 
basic  vortex/dipole  principles. 

Turbulence 

One  of  the  major  differences  between  potential  flow  theory  and  real  fluids  is  turbu¬ 
lence.  Interestingly,  turbulence  by  itself  does  not  affect  significantly  the  time-average 
work  output  of  most  turbomachines.  It  seems  that  the  dynamic  behavior  of  the  turbulent- 
flow  flow  field  is  represented  properly  by  a  Sear’s  function  type  of  potential  flow  force 
and  work  response.  Therefore  this  primarily  linear  three-dimension^  behavior  is  space 
and  time  averaged  by  the  blades  to  give  the  same  overaU  energy  n-ansfer  as  with  steady 
flow.  For  the  acoustic  output,  the  answer  is  different.  As  demonstrated  earlier  in  Eq.  (42), 
the  far-field  acoustic  wave  depends  on  the  time  derivative  of  flow  field  parameters.  Thus 
because  the  nature  of  turbulence  is  unsteady  flow,  the  acoustic  output  from  a  turbulent 
flow  is  fundamentally  different  from  that  of  a  steady  flow.  To  illustrate,  for  a  blade  row, 
stationary  or  moving,  the  time  variations  associated  with  inflow  turbulence  produce  lift 
fluctuations  that  radiate  as  uncorrelated  impulses.  These  have  only  a  limited  similarity  to 
the  deterministic  flows  analyzed  previously. 

Stationary  Blades 

A  second  significant  difference  from  potential  flow  is  that  rotationality  or  distrib¬ 
uted  vorticity,  is  often  present  in  the  fluid,  especially  when  viscosity  and  boundary  layers 
are  important.  For  example,  if  the  inflow  to  a  blade  row  contains  vorticity  normal  to  the 
mean  flow,  the  turning  action  of  a  blade  row  can,  by  conservation  of  angular  momentum, 
convert  this  into  secondary  vortex  motion  with  axis  in  the  flow  direction,  as  shown  in  Fig. 
17.  This  is  of  interest  since  the  inflow  to  a  cascade  usually  contains  boundary  layers  on 
the  inner  and  outer  walls,  with  oppositely  directed  vorticity.  Then  two  oppositely  rotating 
axial  vortices  will  form  at  the  inner  and  outer  ends  of  the  blade  span.  These  vortices  will 
then  be  convected  downstream  with  the  mean  flow.  In  addition,  the  inflow  boundary  layer 
velocity  profile  approaching  a  right-angle  blade- boundary  wall  comer  junction  will  cause 
the  formation  of  a  horseshoe  vortex  as  the  flow  splits  to  pass  on  either  side  of  the  blade. 
Further,  if  the  bound  vortex  distribution  is  non-u^orm  ^ong  the  blade  span,  vorticity 
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will  be  shed  irom  the  trailing  edges  and  a  sheared  flow  will  exist  downstream.  Finally,  if 
the  rate  of  pressure  rise  or  the  turning  through  a  blade  row  is  too  great  in  the  flow  direc¬ 
tion.  the  newly-formed  blade  boundary  layer  can  experience  reverse  flow  and  large 
separated-flow  vortex  strucmres  can  be  generated  in  an  unsteady  fashion.  All  of  the 
above  can  occur  for  stationary  and  moving  blade  rows,  as  shown  in  Fig.  18.  The  result  is 
that  the  downstream  flows  of  real  fluids  can  have  considerable  new  vorticity,  which  may 
be  inherently  unstable,  leading  to  losses  and  sources  of  sound.  Moreover,  for  separated 
flows  where  there  may  be  large  dynamic  pressure  forces  appearing  on  the  blades,  strong 
acoustic  dipole  radiation  can  occur  from  each  impulsive  event. 


Fig.  17.  Conversion  of  inlet  vorticity  by  a  turning  vane  row. 


Rotating  Blades 

When  blades  rotate,  additional  secondary  flows  appear  with  real  fluids  because  of 
the  unique  action  of  the  moving  boundary  with  viscous  surface  layers.  First,  the  boundary 
layer  on  a  rotating  blade  is  affected  more  by  the  centrifugal  force  than  the  through-flow. 
Depending  on  the  machine  geometry  and  surface  pressures,  this  force  field  tends  to  set  in 
motion  a  large  channel  vortex  system  with  axis  in  the  mean  flow  direction.  Secondly,  for 
blades  with  significant  tip  clearance,  some  of  the  bound  vortex  system  can  be  shed 
around  the  tip,  resulting  in  a  region  of  high  shear  at  the  wall.  (For  the  ideal  internal  flow 
machine,  where  the  blade  is  loaded  uniformly,  this  vortex  is  shed  mathematically  into  the 
outer  boundary  and  does  not  cause  complex  flows.)  Also  important  for  the  moving  blade 
case  is  the  fact  that  the  boundary  layer  leaves  the  trailing  edge  in  a  manner  different  than 
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for  a  stationary  blade,  a  fact  that  is  often  overlooked.  This  is  because  when  viewed  in  the 
absolute  frame,  as  in  Fig.  19,  the  boundary  layer  fluid  is  attached  to  a  moving  surface  and 
represents  a  high  energy  flow  moving  at  some  angle  to  the  bulk  fluid.  Clearly  the  details 
of  the  trailing  edge  dynamics  in  this  case  are  different  than  for  flow  over  a  stationary 
blade.  Finally,  it  has  been  observed  that  for  certain  flow  conditions  turbomachines  devel¬ 
op  large  vortex  structures  outside  of  the  blade  row  that  rotate  in  the  same  direction  as  the 
shaft,  but  at  some  lower  speed.  These  vortex  systems  are  called  “rotating  stall”  and  can 
produce  very  strong  dynamic  effects  at  the  inlet  or  the  exit. 


Fig.  18.  Combination  of  blade  row  vortex  systems. 

All  of  the  above  effects  represent  a  transfer  of  fluid  energy  into  new  vortex  or  shear 
layer  flows  from  which  that  energy  cannot  be  recovered  in  any  practical  manner  More¬ 
over,  these  complex  flows  tend  to  be  naturally  unstable  and  to  be  dispersed  by  turbulence, 
evenhially  being  dissipated  by  viscosity.  The  literature  on  turbomachinery  contains  exten¬ 
sive  studies,  experimental  and  theoretical,  on  these  “secondary”  flows  since  they  are  so 
important  to  machine  efficiency.  As  observed  by  Atassi  in  a  recent  review  of  the  Sear’s 
function  properties, 

“The  circulation  aroimd  the  airfoil  changes  in  response  to 
every  unsteadiness  in  the  flow.  For  every  change  in  cir¬ 
culation,  a  vortex  must  be  shed  at  the  trailing  edge  of  the 
airfoil  and  is  carried  away  by  the  mean  flow.  Therefore 
vortices  shed  in  the  wake  represent  a  recorded  history  of 
the  airfoil  non-uniform  motion.” 
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Rg.  19.  Trailing  edge  absolute  velocities  tor  stationary  and  moving  altloiis. 

Although  stated  for  unsteady  flow,  this  observation  applies  to  all  the  shed  vortex 
systems,  since  each  produces  its  own  wake  signature.  It  is  just  this  effect  that  the  Lighthill 
vorticity  moment  integral,  Eq.  (28),  is  intended  to  capture  for  a  single  blade  surface.  For  a 
complete  turbomachine,  which  on  the  average  introduces  torque  and  thrust,  the  result  in 
the  wake  is  a  history  of  impulses  from  each  blade  and  vane,  along  with  any  duct  bound¬ 
ary  conditions,  to  yield  a  total  flow  field  with  changed  energy  and  momentum. 


VORTICES  AND  ACOUSTICS 

From  the  acoustic  standpoint,  the  real  turbomachine  generates  primarily  the  work- 
related  vortex  system  predicted  from  potential  flow  theory.  But  in  addition,  the  set  of 
secondary  vortices  discussed  above  appears  in  a  complex  marmer  that  is  difficult  to  cate¬ 
gorize  in  terms  of  total  fluid  acoustic  behavior.  One  approach  to  assessing  the  overall 
imsteady  acoustics  is  to  assume  that  each  blade  can  be  described  analytically  by  a  time 
and  space  distribution  of  externally  applied  volume  and  force  sources.  Thus  the  mass  and 
momentum  conservation  equations  become: 

^  +  V  ■QV  =  m{x,t)  (100) 

dr 

+  (Vghv  =  Vp+M^^V+/(x,t)  (101) 
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For  flows  where  the  density  changes  and  Mach  numbers  are  small,  the  time  derivative  of 
Eq.  (100)  and  the  divergence  of  Eq.  (101 )  can  be  combined  to  give; 

=  I  WV2)+f  X  +  V  /(102) 

(r  dr  \  Qo  J  dt 


In  this  inhomogeneous  form  of  the  wave  equation,  the  terms  on  the  right  represent  dis¬ 
tributed  acoustic  sources.  If  the  source  activity  is  limited  to  a  small  region  near  the  origin, 
and  the  intervening  fluid  is  uniform  and  stationary,  then  a  Helmholtz  integration,  with  the 
delay  times  retained,  over  the  active  fluid  region,  and  assuming  acoustically  compact  sur¬ 
faces,  gives 


Pa 


1 

4jtro 


Qo^y  •  (V(VV2)  +T  X  V-  v/p„V2\0]^r/c^+  (103) 


Equation  (103)  is  a  general  result  for  far-fleld  radiation  from  a  small  source,  but  to 
appreciate  its  meaning  it  is  useful  to  consider  some  simple  cases. 

Rigid  Sphere 

First,  let  the  source  be  the  rigid  sphere  oscillating  in  an  inviscid  fluid.  The  first  vol¬ 
ume  integral  is  to  be  taken  over  all  the  fluid  outside  the  sphere  using  a  singly-coimected 
region  that  contains  no  singularities.  Doing  this  for  the  case  of  low  Mach  number  results 
in  the  integral  being  negligible  for  an  irrotational  and  inviscid  flow.  The  second  integral 
involving  the  divergence  of  the  force  vector  can  be  converted  to  a  propagating  time  deriv¬ 
ative  by  recognizing  that  in  the  far-fleld; 

(104) 

^  c  df 

Taking  the  second  integral  over  the  surface  of  the  sphere  gives  the  net  force  vector  with 
the  correct  time  delay: 

j  I  fi>-r/c  ^  •  j  j \P^sl-rJc  ds  =  ~  K  •  (105) 

Examination  of  the  last  integral  of  Eq.  (103)  shows  that  it  represents  the  mass  intro¬ 
duction  to  the  region,  which  for  a  constant  volume  sphere  will  be  zero.  However,  as 
discussed  earlier  the  linear  motion  of  the  sphere  is  simulated  by  mass  introduced  at  the 
front  and  removed  at  an  equal  rate  at  the  back.  The  acoustic  time  delay  operation  on  this 
behavior  results  in  the  time  derivative  of  the  first  moment  of  the  surface  integral  of  mass 
addition  in  the  direction  of  the  observer. 
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Therefore  the  total  far-field  acoustic  pressure  for  a  moving  solid  sphere  is; 


-  -  "no-  d  =r  dU 

=  - —  FA  +  gJi— 

AjtToC  dt  dt  f 


AjtToC  dt 


This  can  be  recognized  as  the  combined  added  and  displaced  mass  terms  in  the  Lighthill 
equation  for  the  force  on  a  moving  rigid  body.  (Eqs.  (26a)  and  (26b)  or  Eq.  (48)).  Inter¬ 
estingly,  these  two  force  vectors  need  not  be  in  the  same  direction,  since  the  added  mass 
term  responds  to  the  relative  velocity. 

Moving  Airfoil 

Taking  next  the  case  of  a  moving  airfoil  of  zero  thickness,  the  volume-related  terms 
are  neglected,  but  the  integration  of  Eq.  (103)  must  now  account  for  the  divergence  of  the 
work  force  vector  due  to  the  relative  flow  over  the  airfoil.  This  can  be  computed  as  above 
so  that: 

f  [  f  [V,  f^-l-r/cdV =•^^11  * = T  ■  t?-W/I08) 

For  a  moving  airfoil,  the  vorticity  in  the  first  integral  of  Eq.  (103)  is  somewhat  more  dif¬ 
ficult  to  interpret  since  the  work  force  can  produce  vorticity  in  the  fluid.  When  there  is 
vorticity  distributed  in  the  fluid,  it  may  be  recalled  that; 

f  {^dA^Tny  (27b) 


where  the  surface  integral  is  taken  such  that  the  vorticity  is  normal  to  the  surface,  and  the 
singularity  strength  T  represents  all  the  vorticity  in  that  area.  If  now  the  local  shed  vortic¬ 
ity  vector  is  aligned  with  the  local  velocity,  as  for  a  fixed  blade,  or  if  there  is  no  shed 
vorticity  in  the  fluid,  as  for  a  blade  with  spanwise  uniform  circulation,  the  first  volume 
integral  is  again  negligible.  Alternately,  if  the  blade  is  of  finite  span,  and  moving,  the 
shed  vorticity  will  include  both  positively  and  negatively  rotating  parts.  These  will  be  at 
some  angle  to  the  velocity  and  must  be  present  to  account  for  the  span-wise  enthalpy  dis¬ 
tribution.  As  a  consequence,  over  the  volume  of  integration  both  positive  and  negative 
span-wise  force  components  will  appear  in  the  wake.  Once  introduced  into  an  ideal  fluid, 
this  pattern  of  opposing  forces  will  not  change,  except  for  convection.  Therefore  it  will 
not  produce  a  significant  first-order  time  derivative  when  Eq.  (104)  is  applied  to  the  vol- 
tune  integral  of  the  vorticity-velocity  cross  product.  Hence  for  an  airfoil  of  negligible 
volume,  evaluation  of  Eq.  (103)  yields  only: 


where  the  Lighthill  vorticity  moment  integral  is  to  be  taken  over  all  the  volume,  and  in¬ 
cludes  the  internal  bound  and  the  external  shed  vorticities. 

Combining  the  above  effects  for  a  basic  ideal  turbomachine,  the  total  dipole  acous¬ 
tic  sources  associated  with  one  rigid  blade  are  found  to  be  due  solely  to  forces  applied  to 
the  fluid  from  without. 
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^(0=1  \ ~^Yt\  1 (110) 


These  forces  cause  the  expanding  vortex  loops  necessary  when  there  is  work  addition  to 
the  fluid  and  include  the  primary  and  secondary  flows.  They  also  represent  the  motion  of 
the  fluid  mass  displaced  normal  to  the  blade,  along  with  the  blade  surface  vortex  loops 
necessary  to  produce  the  correct  surface  tangential  velocity.  In  typical  turbomachine  de¬ 
signs  the  volume-related  term  will  be  small  when  compared  with  the  work  term. 

Turbulent  Sources 


Qearly,  in  regions  of  the  flow  where  there  are  no  surfaces,  rti  =f-  o,  and  the  only 


acoustic  sources  will  be  the  natural  fluid  instabilities,  combined  with  turbulence.  In  this 
case  only  the  first  integral  of  Eq.  (103)  is  needed  to  represent  the  dynamic  forces  in¬ 
volved  for  internally  interacting  fluid  motions.  The  terms  are  the  vector  form  of  the 
well-known  Lighthill  stress  tensor  representation  of  the  same  phenomena,  i.e.. 


( 


V(VV2)  +  ?  X  V-H 


Qo  j 


d\QoUiUj-^Xij)  ^  d^Tjj 
dyidyj  dyidyj 


(111) 


Thus,  away  from  surfaces  the  instabilities  in  real  fluids  produce  an  assembly  of  opposing 
motions  that  vary  with  time  in  their  magnitude,  location  and  direction.  In  these  motions  it 
is  unclear  how  to  distinguish  between  kinetic  energy  and  vorticity  transport  terms,  hence 
the  utility  of  the  stress  tensor  form  that  accounts  for  both.  If  the  inflow  to  the  machine  is 
naturally  turbulent  then  these  local  velocity  fluctuations  will  cause  fluctuating  pressure 
distributions,  and  an  altered  pattern  of  wake  vortices  as  described  by  Atassi.  In  addition, 
fluctuating  flow  across  the  solid  body  will  cause  added  mass  impulse  pressures  to  appear. 
Figure  20  illustrates  these  acoustic  sources  for  a  blade  row.  It  is  possible  that  these  un¬ 
steady  sources  can  be  clarified  by  CFD  methods  such  as  large  eddy  simulation  or  direct 
numerical  simulation  of  turbulence. 

In  summary,  many  vortex  systems  occur  for  flows  through  turbomachinery,  and 
many  levels  of  analysis  exist  to  understand  and  control  these.  The  combined  techniques 
of  potential  and  real  fluid  analysis,  as  characterized  by  the  many  contributions  of  Howe 
seem  to  be  an  effective  way  to  develop  a  building  block  approach  to  acoustic  perform¬ 
ance.  Suffice  it  to  say  that  the  first  acoustic  goal  should  be  to  avoid  the  generation  or 
propagation  of  tonal  signals  associated  with  shaft  rotation.  The  second  acoustic  require¬ 
ment  should  be  to  understand  the  sources  of  broadband  noise  and  to  develop  iimovative 
designs  that  suppress  these  sounds. 
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Fig.  20.  Surface  force  dipoles  and  bulk  flow  quadrupoles. 


CONCLUSIONS 


The  principles  of  unsteady  potential  flow  theory  provide  a  convenient  and  unified 
approach  to  understanding  the  work  and  acoustic  properties  of  turbomachinery.  The  clas¬ 
sic  mathematics  of  Lamb,  augmented  by  the  contemporary  work  of  Lighthill,  and 
combined  with  specific  applications  by  other  investigators  provides  a  reasonably  com¬ 
plete  analytical  description  of  ideal  turbomachines.  From  this  it  can  be  seen  that  all  of  the 
energy  transfer  functions  of  rotating  machinery  are  related  to  moving  and  expanding 
three-dimensional  closed  loop  vortex  systems  originating  from  the  torque  and  thrust 
forces  on  each  discrete  blade.  Although  the  internal  details  of  the  flow  field  differ  be¬ 
tween  axial  and  centrifugal  machines,  the  forces  can  be  generalized  by  considering  their 
origin  on  each  blade,  and  the  combination  of  their  shed  vortex  loops  to  produce  the  total 
enthalpy  change  in  the  external  flow. 


In  the  mathematical  models  for  these  unsteady  potential  flows  there  exists  a  near- 
field  solution  that  accounts  for  the  work  done,  and  a  far-field  solution  representing  the 
acoustic  behavior.  Because  of  the  mathematical  equivalence  between  vortex  loops  and 
combinations  of  dipoles,  it  is  straight  forward  to  show  how  the  acoustics  are  derived  from 
the  work-related  flow  fields.  The  basic  acoustic  result  for  propagation  into  a  three-dimen¬ 
sional  field  due  to  a  dynamic  force  applied  near  the  origin  is  the  dipole  pressure  given  by: 
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Fortunately,  at  geometrically  remote  points  a  single  force  can  account  for  the  work  and 
volume  displacement  effects  of  a  single  blade  moving  through  a  fluid  field.  That  force 
can  be  computed  from  a  knowledge  of  the  pressure  distribution  on  the  blade,  or  from  an 
expression  developed  by  Lighthill  that  uses  kinematic  flow  variables  from  the  entire  fluid 
field; 

^0  =  I  j Pn4s  =  Y^rl  j  ^ 

Even  though  the  initiating  force  may  be  well-defined,  solutions  for  the  far-field 
acoustics  are  complicated  by  the  need  to  account  properly  for  the  delay  time  operation. 
Because  of  the  circular  symmetry  of  all  rotating  machinery,  ideal  solutions  for  the  acous¬ 
tic  field  introduce  Bessel  functions  as  shown  first  by  Gutin.  For  rotor-stator  combinations 
in  a  free  field  with  equally  spaced  blades  and  vanes  the  Bessel  function  order,  m,  is  deter¬ 
mined  precisely  by  the  following  integer  relations, 

m  =  nB-^kV  axial  propagation 
m  =  nB+kV+ I  radial  propagation 


When  the  Bessel  function  order  is  zero,  for  propagation  in  the  directions  indicated,  strong 
radiation  can  occur  due  to  geometric  alignment  of  the  blade  and  vane  harmonic  patterns. 
The  radiated  sound  is  always  at  blade  passing  frequency  and  its  harmonics,  while  Bessel 
functions  control  the  spatial  pattern. 


For  acoustic  radiation  in  a  duct,  the  well-known  results  of  lyier  and  Sofrin  apply 
with  axial  propagation  occurring  only  when  the  blade  harmonic,  Bessel  function  order 
and  relative  mach  number  satisfy: 

nB 

—  Mret  >  1  Duct  (axial)  propagation 


Even  for  low  speed  flows,  if  rotor  stator  combinations  allow  m  =  o,  those  harmon¬ 
ics  will  radiate  in  the  plane  wave  mode.  If  this  Mach  number  criterion  can  be  met  for  m  ^ 
o,  the  propagating  wave  is  in  the  form  of  a  “spinning  mode”. 

These  ideal  propagation  analyses  predict,  for  a  low  Mach  number  multi-bladed  rotor 
in  a  uniform  flow,  that  tonal  radiation  should  be  very  small.  But,  for  rotor/stator  or  stator/ 
rotor  interactions  the  tonal  acoustics  can  be  quite  strong.  Nevertheless,  such  tones  can  be 
suppressed  by  using  the  well-known  design  nile; 

V  >2B 


which  has  the  effect  of  preventing  low  order  Bessel  functions  with  potentially  large  am¬ 
plitudes  from  occurring  in  the  solution. 

These  precise  mathematical  functions  provide  clear  rules  for  controlling  the  tonal 
acoustics  of  turbomachines.  However,  there  remain  sources  of  sound  that  are  broadband, 
and  yet  to  be  characterized,  but  which  may  be  tmderstood  better  using  vortex/dipole  con¬ 
cepts. 

Extension  of  the  unsteady  vortex/dipole  principles  to  turbomachines  handling  real 
fltiids  with  turbulence  and  secondary  flows  provides  some  additional  useful  insights  into 
the  work  and  acoustics.  For  example,  while  the  time-average  work  function  due  to  blade 
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motion  is  relatively  insensitive  to  inflow  turbulence,  the  acoustic  behavior  is  quite  sensi¬ 
tive  since  it  is  precisely  the  time  derivative  of  the  unsteadiness  that  provides  the  acoustic 
source.  In  addition  blades  in  real  machines  generate  a  number  of  secondary  vortex  sys¬ 
tems  that  are  generally  unstable  and  these  result  in  unwanted  mechanical  losses,  as  well 
as  acoustic  signals  associated  with  the  time  derivatives  of  their  growth  and  decay.  Un¬ 
steadiness  in  the  secondary  flow  that  occurs  away  from  the  influence  of  blade  surfaces 
will  produce  quadrupole  acoustic  sources  as  characterized  by  the  Lighthill  stress  tensor. 

The  applications  of  unsteady  potential  flow  theory  presented  herein  are  believed  to 
represent  a  unified  and  consistent  description  of  the  near  and  far-field  effects  for  low 
speed  turbomachines.  Each  topic  was  selected  to  illustrate  to  a  particular  fundamental 
point  and  to  show  how  that  understanding  contributed  to  the  larger  view.  As  stated  initial¬ 
ly,  these  results  are  focused  on  cases  where  the  stationary  and  moving  surfaces  are  rigid, 
and  where  the  acoustic  source  is  compact.  It  is  thought  that  this  introduction  will  provide 
a  sound  basis  for  future  machine  development. 
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